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SUMMARY 


This study focused on characterizing matrix cracking and delamination 
behavior in multidirectional laminates. Static tension and tension-tension 
fatigue tests were conducted on two different layups. Damage onset, 
accumulation, and residual properties were measured. Matrix cracking was shown 
to have a considerable influence on residual stiffness of glass epoxy laminates, 
and could be predicted reasonably well for cracks in 90° plies using a simple 
shear lag analysis. A fracture mechanics analysis for the strain energy release 
rate associated with 90° ply matrix crack formation was developed and was shown 
to correlate the onset of 90° ply cracks in different laminates. The linear 
degradation of laminate modulus with delamination area, previously observed for 
graphite epoxy laminates, was predicted for glass epoxy laminates using a simple 
rule of mixtures analysis. The strain energy release rate associated with edge 
delamination formation under static and cyclic loading was difficult to analyze 
because of the presence of several contemporary damage phenomena. 
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Surface area of a matrix crack 
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Multiple linear regression parameters 

0/90 interface delaminated fraction 

+45/-45 interface delaminated fraction 

Laminate half-width 

Matrix crack density 

90° ply matrix crack density 

45° ply matrix crack density 
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Sublaminate modulus 
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Strain energy release rate 

Critical strain energy release rate for damage onset 
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INTRODUCTION 


The application of composite materials in modern helicopters is being 
extended to more complex and safety-critical components which require a deeper 
understanding of the material behavior to guarantee reliability. In the future, 
civil and military authorities may require a failsafe or damage-tolerant design 
of composite helicopters. This need is recognized in principle by the helicopter 
manuf acturers , but it is clear that material properties, test methods, and 
application of composites to real components must be further investigated in 
order to achieve reliable structures. Because glass/epoxy composites are already 
being used extensively in helicopter rotor systems, an investigation of the 
fatigue behavior of glass/epoxy laminates was undertaken in the hopes of laying 
a foundation for future damage assessment of helicopter structural components. 

Composite laminates commonly employed in helicopters, as in many other 
applications, usually show evidence of damage before failure. The usually 
adopted "first ply failure" criterion [1] can very roughly predict transverse 
ply matrix cracking but cannot predict delamination. Both of these two damage 
modes may occur early in the life of the component but may have only a small 
effect on overall stiffness and strength properties. 

The scope of this work is to study these two phenomena from a fracture 
mechanics point of view, offering, as far as possible, some simple but reliable 
criteria to predict onset and accumulation of matrix cracking and delamination. 
Extensive work has already been done with good results using fracture mechanics 
to model both delamination [ 2-5] and matrix cracking [5-13], but only in a few 
cases are these results readily applicable to design. 
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In this work a new approach is developed for characterizing matrix cracking 

onset and growth. The critical strain energy release rate for formation and 

accumulation of matrix cracking in the 90° plies of [0/90] and [±45/0/90] 

s s 

glass/epoxy composites were calculated and plotted for a range of crack 
densities in the transverse plies. 

A previously developed closed form model is available to predict edge 
delamination behavior in composite laminates [2]. In the present study this 
model was applied to glass/epoxy quasi-isotropic laminates, where the presence 
of pre-existing matrix cracking was observed. The theoretical model was then 
compared with experimental results. Critical strain energy release rates 
associated with edge delamination were measured for the material studied. 

EXPERIMENTS 

Two E glass/epoxy laminates were tested in quasi-static tension and in 
tension-tension fatigue. The material was T1 90/X751 /50^ 1 ^ with an average ply 
thickness of .225 mm, 50$ fiber volume, 125 °C curing temperature. The measured 
tensile failure strain in the fiber direction was 1 .76$ and the following 
elastic properties were measured from unidirectional specimens: 

E 1 = 43.7 GPa E 2 = 16.0 GPa 

G 12 = 6.3 GPa v 12 = .284 

d ^Produced by Cyanamid Fothergill Ltd., U.K. 
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The layups tested were [0/90] and [±45/0/90] . The first sequence was 

s s 

chosen to isolate the 90° ply matrix cracking from other damage phenomena; the 
second sequence was chosen to highlight the edge delamination at the 0/90 
interface, because of the high peeling stresses which are generated [2]. Tests 
were run on an Instron servo-hydraulic machine. Axial strain was derived from 
displacements measured by a pair of LVDT's over an 82 mm gage length. A high 
intensity light was shone through the specimen to observe and photograph the 
damage. Matrix crack density and delamination were measured from high contrast 
photographic prints like the one shown in fig. 1. For each laminate five quasi- 
static tension tests were run with a ram speed of 0.5 mm/min. At specified loads 
the specimens were unloaded to record damage and residual stiffness (tangent 
modulus). Tension-tension load-controlled fatigue tests were conducted with a 
repeated sinusoidal load of 6 Hz frequency and 0.1 load ratio; 20 fatigue tests 

for the [0/90] lay-up and 25 tests for the [+45/0/90] lay-up were performed, 
s s 

covering the whole range of possible fatigue lives from 1 to 10 6 cycles. The 
cyclic loading was interrupted at specified intervals, to monitor damage and 
measure static stiffness. 

DAMAGE OBSERVATIONS 


Figure 2 shows schematically the damage types observed on quas i- i sotrop ic 
specimens, namely 90° ply matrix cracking, edge delamination at the 0/90 
interface and local delaminations originating from the angle ply cracks. 
Depending on the type of test ( static or fatigue ) and on the load levels, 
damage accumulation showed some qualitative and quantitative differences. 

From previous works it is known that matrix cracks do not grow indefinitely 
in number, instead they reach a "saturation value", which seems to be a function 
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of the laminate characteristics [8J. In static tests failure occurred when the 


number of 90° ply matrix cracks was still low and angle-ply cracks had not yet 
developed. The edge delamination at failure was extended over 20% of the 
interfacial area and the stiffness loss was typically about 15$ of the initial 
stiffness. 

Fatigue damage accumulation is shown in fig. 3. This figure shows the 
results of a single test run at 1 60 MPa maximum sinusoidal load (R=0.1), which 
is about 60$ of the laminate static failure stress, and is typical for damage 
developed during fatigue tests, run on quasi isotropic glass epoxy laminates . 

The build up of the laminate damage is plotted against fatigue cycles. 
Figure 3(a) shows the matrix cracking formation in 90°, +45° (external) and -45° 
(internal) plies. For each of these plies, matrix cracks reach a saturation 
level at different number of cycles. The cracks in the 90° plies accumulate 
faster than matrix cracks in the others, so that all other damage may be 
considered to happen after the 90° plies damage reached saturation. Fig. 3(b) 
displays the edge delamination and local delaminations growth. Comparing 
fig. 3(a) and 3(b), the delaminations and cracking of +45 and -45 plies appear to 
be almost contemporary, so that their effects must be superimposed. The 
cumulative effect of the damage on the specimen longitudinal modulus E/E 0 is 
shown in fig 3(c). 

The same type of tests were performed on [0/90] laminates in order to 

s 

isolate the behavior of 90° ply matrix cracking. The [0/90] laminate showed 

s 

only 90° matrix cracking in static tests. Also in fatigue, the damage was mainly 
90° matrix cracking, but the onset of matrix cracking in 0° plies, i.e. 
longitudinal splits, was observed before final failure. Total stiffness loss 
just before failure was 15$ in static tests and more than 20$ in fatigue tests. 
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ANALYSIS OF MATRIX CRACKING 


Matrix cracking occurs when tensile stresses normal to the fibers are 
present. On glass reinforced composite materials this kind of damage is 
responsible for a considerable fraction of the total stiffness loss both in 
static and in fatigue tests. To predict the matrix cracking effect on different 
laminates and loadings, available models have been studied and compared to 
experimental results, and a fracture mechanics approach has been applied to 
characterize matrix crack onset and growth. 

Stiffness loss in a [0 /90 ] laminate 

m ns 

A closed form expression for stiffness loss in a [0 /90 ] laminate as a 

m ns 

function of crack density D (fig. 4) was derived previously from a shear lag 
analysis [13] as 


E 

E 


X 

o 


1 


. + nE, tanh (A/2D) 
mE, X/2D 


( 1 ) 


where 

X 2 . 3 G 12 E Q (n+m) 

h 7 E t E 2 n 1 IT 

and E^, E 0 are the moduli of the cracked and the uncracked laminate 
respectively. The E lf E 2 and G l2 are the lamina longitudinal, transverse and 
shear moduli respectively, and h is the ply thickness. 
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To evaluate the degree of approximation introduced by eq.(1) in the 
laminates presently studied, this equation has been compared with results from a 
NASTRAN finite element model (see Appendix A) and with experimental data. Fig. 5 
compares equation (1) (continuous line) and Finite Element Model (diamonds and 
dashed line) with data obtained in static and fatigue tests on [0/90] 
laminates. The small difference observed between the finite element results and 
eq.(1) will produce consequences that will be discussed later. The comparison 
with the experimental data also confirmed the validity of equation (1) for 
elevated crack densities. 

Stiffness Loss in Other Laminates 

Although equation (1) is only valid for [0/90^ stacking sequences, its 

use can be extended to model stiffness loss due to matrix cracking in the 90° 

plies, surrounded by 0° laminae, in other laminates of arbitrary orientation. 

Previously, a technique was developed for estimating the reduction in for a 

cracked 90° ply as a function of the crack density [6]. This technique is 

depicted in fig. 6. First, the reduction in laminate modulus, E x /Eq, of a [0/90] s 

laminate is determined as a function of the crack density in the 90° ply (fig. 

6(a)). In the current study, this was determined using eq.1 (fig. 5). Next, the 

linear relationship between the reduced value of E 2 in the cracked 90° ply and 

the [0/90] laminate modulus is calculated from laminated plate theory (fig. 
s 

6(b)). The reduced laminate modulus, E x /E Q , calculated when the effective 
transverse lamina modulus in the 90° ply, E^/E^, is set to zero is designated 
E # /E q , and is the value obtained from the ply discount method where the cracked 
90° ply is assumed to carry no load. Finally, the equation of the straight line, 
relating E£/E 2 to E x /E Q (fig. 6(b)), is combined with the equation describing 
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E x /E q as a function of crack density (fig. 6(a)) to generate the equation for a 

master curve relating the effective transverse modulus of a cracked 90° ply in 

an arbitrary laminate, E£/E 2> to the crack density in that 90° ply (fig 6(c)). 

This master curve may then be used to determine the appropriate reduction in E 2 

for a cracked 90° ply in another laminate, such as the [45/— ^5/0/90j laminate 

s 

tested in the current study. 

In ref. 6 it was assumed that the longitudinal modulus of a laminate with 

matrix cracks is a linear function of the " effective damaged lamina " 

transverse modulus E' normalized to the undamaged lamina transverse modulus E 2 

2 

(fig. 6(b)). The equation of this straight line yields: 


E 2 '(D,w) E (D, w) - E 

= — - ( 2 ) 

E 2 E 0 - E # 

In reference 6, the finite element method was used to determine the [0/90] 

s 

laminate stiffness loss as a function of crack density, E x (D,w), whereas in the 
present study it was possible to start from a closed form expression (eq.1). 
Equation (2) may be rearranged to yield 

E x • h * <E » - V ( 3> 

A precise expression of eq. 3 may be obtained through laminated plate 

theory (App. B), but the result shows no significant difference from eq. 3* 

Eq. 2 gives E'/E 2 as a function of crack density and ply thickness, w=2nh 
2 

(fig. 4), for a laminate stacking sequence where the relationship between E x and 

D*w is known. In the present case, a [0/90] glass/epoxy laminate was modeled 

s 
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using equation (1). This approach is valid only if the 90° ply is constrained by 

adjacent 0° plies. For other cases, where more complex shear stresses are 

involved, the relationship between E^ and crack density may be obtained from 

finite element analysis and substituted into eq.2 to determine the functional 

relationship between E^/E^ and crack density. 

To check the validity of representing the effect of transverse ply cracking 

on laminate stiffness in terms of a calculated effective modulus E' 2 (D,w), the 

longitudinal modulus of a [±45/0/90] laminate was calculated and plotted in 

s 

fig. 7 (solid line) for several 90° ply crack densities using E^(D,w) for the 
cracked 90° plies obtained from eq.2. The quasi-isotropic laminate containing 
various crack densities was also modeled using the NASTRAN finite element code 
(App. A). The laminate stiffnesses calculated by FEM are shown in fig. 7 as 
diamond symbols connected by a dashed line. The agreement between the calculated 
stiffness loss, from eq.2 and FEM, and the experimental data (triangles in fig. 
7) was fairly good. 
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Strain Energy Release Rate 


A fracture mechanics analysis was applied to develop a tool capable of 
predicting onset and growth of matrix cracking. For an elastic body containing a 
flaw of area A that extends under a constant nominal strain, e, the strain 
energy release rate, G, may be written as 


G 


-V 


e 2 dE 
2 dA 


(4) 


where V is the volume of the body and dE/dA is the rate of stiffness change as 
the flaw extends [2]. 

The fundamental assumption of this work is that all of the N matrix cracks 
in the volume V may be treated as a single equivalent flaw of area A, that is 


A = I 


N 

i=1 


A. 

l 


(5) 


and that all the matrix cracks have the same length, w, equal to the thickness 
of the cracked laminae (fig. 4). Therefore, for a laminate of width, 2b, length, 
%, and total thickness, t, 


A = N 2b w, and V = 2b % t 


( 6 ) 


Because N = i D, A can be espressed in terms of crack density D: 


A = 2b H w D, and dA = 2b l w dD 
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With these assumptions G may be written as 


t e 2 dE 
w 2 dD 


( 7 ) 


Differentiating eq. 1 and substituting into eq. 7 yields a closed form 
expression of the strain energy release rate for matrix cracking in a [0^/90^^ 
laminate: 


tanh( 1/2D) 1 

A/2D cosh 2 ( A/2D) 

(n+m) E ? E fl e 2 r ] (8) 

m E, D 2 L J 

ri + nE ? tanh(A/2D) -, 2 
U + mEj A/2D J 

The closed form expression is particularly suitable for comparing different 

lay-ups and materials. Fig. 8 shows the calculated effect on strain energy 

release rates of different lay-ups and reinforcing fibers. In the case of Glass 

fibers (solid and dashed curves), the effects of doubling either the 90° or the 

0° ply thickness are shown. Doubling the 90° ply thickness results in a doubling 

of the strain energy release rate, as it is evident comparing the curves for the 

[0/90] and [0/90 o ] laminates. Doubling the number of 0° plies, however, has 
s 2 s 

very little effect on G as seen when comparing the curves for the [0/90^^ and 

[0 o /90 o ] laminates. The difference in the G curves for these two laminates is 
2 2 s 

greatest at small values of crack density. Hence, the onset of 90° matrix cracks 
may occur differently for the two laminates, but the accumulation of cracks 
should be similar for both. Figure 8 also shows a comparison of (0/90) 

s 

glass/epoxy and graphite/epoxy laminates. As expected, graphite/epoxy shows a 
lower and flatter G curve, because 90° ply matrix cracking has a lesser effect 
on laminate stiffness in graphite/epoxy compared to glass/epoxy [2]. 
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Figure 9 compares strain energy release rate as a function of crack 
density, for a [0/90] glass/epoxy laminate, calculated both from eq. 8, and 
from eq. 7 using dE/dD derived from finite element results. The difference in 
the G values from the two techniques results from the difference in slope of the 
stiffness vs. crack density expressions, which can be observed in figure 5. At 
the moment, it is not possible to establish which of the two models is closer to 
the actual behavior in the initial region of low crack density, since both rely 
on some simplifying assumptions, and the scatter in the experimental data shown 
in fig. 5 makes it difficult to resolve which analysis is more appropriate. In 
the high crack density region, however, the finite element analysis seems to 
comply better with the experimental data and with the assumption that for high 
crack density G must reach zero in order to explain the observed phenomenon of 
crack density saturation. 


Matrix crack onset 


To evaluate G for the onset of 90° ply matrix cracking in a [0/90] 
o m n 3 

laminate, it is easier to use a simpler form of eq. 8 with D set to zero 


= (n+m) E, E n 2 
'c m Ej X e c 


(9) 


A similar expression was derived from a shear lag analysis in ref. [10] and [11]. 
Equation 9 differs only in that the factor 3 included in the X expression ( eq . 1 ) 
is replaced by a 4 in ref. [10] and [11]. This difference is a result of the 
different distributions of the longitudinal displacements in the transverse ply 
used in the shear lag analysis. This distribution was parabolic in the 
derivation of eq. 1 [13] and linear in ref. 10. 
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For a generic laminate with 90° plies surrounded by 0° plies, dE/dD is 
easily calculated from eq.(2) when D=0. This calculation introduces an 
approximate linearized expression for effective damaged lamina transverse 
modulus, obtained from eq.(1), as D approaches zero. 

Crack growth 

As crack density increases, eq. 8 gives a lower value of G for a fixed 
strain. In a constant strain fatigue test no more cracks should appear, as soon 
as G is reduced below its critical value. In the present study, where all the 
tests were run under load control, an experimental proof was not possible and, 
as already stated, in most cases the matrix cracks reached a saturation density 
value. The presence of this saturation density is not totally explained with the 
closed form expression, because in eq.8, the dE/dA term never goes to zero. If 
dE/dA does not go to zero, then an increased stress will always generate more 
cracks. The finite elements analysis however, shows more clearly that dE/dA is 
very close to zero when the crack density is near to the saturation value, no 
matter what the stress level (fig. 9). 

Looking at the problem from another point of view, the "apparent" critical 
G value, G D , would increase with increasing crack density. Therefore, a crack 
resistance curve concept was borrowed from classical fracture mechanics [17], 
showing the progressive laminate resistance to formation of new matrix cracks as 
the number of these cracks increases. This increased crack growth resistance is 
not attributed to other energy absorbing damage mechanisms, but results only 
from the constraint and load transfer capability of the neighboring plies 
adjacent to the cracked ply. 
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Crack growth resistance curves (R-curves) for 90° ply matrix cracking were 

obtained from quasi-static tension tests. Five tests of this kind were performed 

on [0/90] and [±45/0/90] laminates. Figure 10 shows G„ calculated from strain 
s s R 

and crack density data for both laminates, using equation (8). The first point 

plotted for each curve corresponds to the onset of the first matrix crack, and 

shows G as would be calculated from eq. 9. The large scatter observed in quasi- 
c 

isotropic laminate data for high values of crack density is mainly due to the 
superposition of other damage modes (i.e. delamination) that can influence the 
laminate behavior. Results from the two laminates are superimposed in fig. 10 
and show good agreement, thus indicating that the R-curve approach is suitable 
for predicting matrix crack growth in 90° plies of a laminate, when the cracked 
90° plies are surrounded by 0° plies of the same material. 


FREE EDGE DELAMINATION 


A model for stiffness loss and strain energy release rate for free edge 
delamination onset and growth in graphite/epoxy laminates has been presented by 
O'Brien [2]. In this model, stiffness decreases linearly as delamination size 
increases and can be expressed by 


E 

x 




( 10 ) 


where E^ is the laminate modulus before delamination, a/b is the ratio of 

* 

delamination size to specimen width, and E is the modulus of the laminate 


completely delaminated along one or more interfaces. 
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The delaminated modulus, E , can be calculated as the sum of the products 
of the moduli and thicknesses of the sublaminates which are created by the 
delamination. That is 


* 

E 



E. 

l 


T 


t. 

1 


( 11 ) 


where t^ and E i are the thickness and the moduli of the i sublaminate and M is 
the total number of sublaminates. The strain energy release rate for edge 
delamination can be evaluated by differentiating eq. 10 with respect to a and 
substituting into eq. 4 as shown in ref. 2, 


G = 


e 2 t 
2 



* 


E ) 


( 12 ) 


The critical value of strain energy release rate G c for edge delamination is 
obtained from the actual onset strain e c - This approach has been demonstrated to 
be useful for evaluating and predicting free edge delamination onset and growth 
in graphite/epoxy laminates [14]. 

For glass/epoxy laminates, the same approach has been attempted. However, 
several differences were noted. In graphite/epoxy, stiffness loss due to 90° ply 
matrix cracking was small [2,6], while this is not the case for glass/epoxy 
laminates [ 9 - 11 , 131 . In addition, because G for 90° ply matrix cracking is 
higher for glass/epoxy than graphite/epoxy laminates with the same layup 
(fig. 8), and because the glass/epoxy has a lower Poisson mismatch, included in 
the (E q - E ) term of eq. 12, and hence a lower driving force for edge 
delamination than graphite epoxy, the glass/epoxy laminate is more likley to 
have extensive matrix cracking before the onset of delamination. These matrix 
cracks influence the stress and strain field in the region where the 
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delamination will grow and they act as a delamination resistance mechanism [18- 
203. Also, in glass/epoxy, because of the lower Poisson ratio mismatch, edge 
delaminations grow only a little in quasi-static conditions before the final 0° 
fiber failure. Furthermore, for the glass-epoxy laminates tested in fatigue, the 
effect of delamination on laminate stiffness loss is partially obscured by the 
simultaneous growth of ±45° matrix cracking, which was also observed previously 
for graphite epoxy [20]. 

As shown earlier in fig. 3, there are many contributions to stiffness loss 
in a [±45/0/90] s glass/epoxy specimen during a fatigue test. To evaluate the 
contribution of each damage mode, a multiple linear regression analysis was 
performed on stiffness loss data. Because the relationship between stiffness 
loss and edge delamination size is linear (eq.10), and because a linear 
approximation between stiffness loss and 90° ply matrix crack density seems 
reasonable (fig. 5), a linear relationship was assumed between stiffness loss and 
angle ply matrix crack density, and 0/90 edge, and +45/-45 local interface 
delamination sizes, as follows: 

a o ai ^90 ~ 3z ^45 ~ a3 ^- 4 c 5 ~ ~ (13) 

The set of coefficients from a 0 to a 5 , for which equation (13) best 
approximated the experimental data, was evaluated statistically. The result of 
this analysis is shown in Table 1A. As previously shown in fig. 3, the 0/90 
interface edge delamination and the angle ply cracking occur simultaneously. 
Therefore, the a 2 , a 3 , and a 5 coefficients are given by the regression analysis 
with a high degree of uncertainty. The contribution of each damage mode to 
stiffness loss is represented by the product of the regression coefficient in 
Table 1A and the appropriate damage parameter in eq. 13 . Therefore, the 
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regression parameters alone do not indicate the relative contribution of each 
damage mode to stiffness loss. For example, the analysis indicated an a 2 value 
that was nearly twice a 3 , but fig. 3(a) shows that the crack density in the -45° 
plies is nearly three times the density in the +45° plies. Hence, the external 
45° matrix cracks had less influence on the laminate stiffness loss than the 
internal -45° matrix cracks. Unfortunately, there is currently no analytical 
model that directly relates stiffness loss to angle ply matrix crack density to 
verify this observation. 

In order to achieve a better experimental evaluation of the edge 
delamination contribution to stiffness loss in the quasi-isotropic laminate, a 
similar statistical analysis was applied to the static test data, where only 90° 
ply cracks and edge delaminations were found. For this special case, eq.13 
reduces to: 


E_ 

E 0 


a 0 



(14) 


The set of coefficients obtained for this three parameter regression analysis is 
reported in Table IB. Dividing eq.10 by E 0 and comparing to eq.14 indicates that 
a 0 = 1 , and that 

a 5 = (E 0 - E*)/E 0 (15) 


Therefore, a direct comparison may be made between a 5 calculated from the 

regression analysis and the value calculated from laminated plate theory using 

eq.15. According to the regression analysis, the delamination contribution to 

* 

the longitudinal modulus loss gives an (E 0 -E ) value of 2500 MPa. However, using 

* 

laminated plate theory to evaluate the terms in eq.15 yields an (E 0 -E ) value 
ranging from 1020 MPa to 770 MPa, depending whether or not the preexisting 90° 
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ply matrix cracks are considered. Hence, G determined using (E 0 -E*) from 
experiments, where significant 90° ply matrix cracks exist ahead of the 
delamination front, would be larger than G determined from calculated 
properties. This observation is consistent with previous studies that showed an 
edge delamination resistance curve could be generated and correlated with the 
increase in 90° ply matrix crack density [2,19,20]. 

These results suggest that the interlaminar stress field near the 
delamination front in the 0/90 interface is influenced by the presence of 90° 
ply matrix cracking. A three dimensional model is required to accurately account 
for the interaction between matrix cracks and delamination. 

FRACTURE TOUGHNESS EVALUATION 

From the analysis above, critical strain energy release rate values were 
obtained for both matrix cracking, eq.(9), and edge delamination, eq.(12), in 
static tension and tension-tension fatigue, and are shown in Table 2. In both 
cases dE/dA values, used to calculate strain energy release rate, were taken 
from experimental results. 

The G q for matrix cracking and edge delamination both decreased 

dramatically in fatigue compared to the static values. This is consistent with 

similar data generated on graphite reinforced composites with a variety of 

matrices [21-23]. The G for delamination was found to be higher than the G for 

c c 

matrix cracking. Although different fracture modes are involved for these two 
mechanisms, this mode difference is not sufficient to explain the large observed 
difference in G c . Other sources [15] report static G Ic for matrix cracking of 
reinforced epoxy composites similar to the measured value reported in Table 2. 
The problem of evaluating G c for delamination onset is again closely connected 
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to the evaluation of the stiffness loss due to delamination and is, therefore, 


affected by the uncertainties reported in the previous section. However a 
comparison with data previously obtained on different composites materials [16] 
showed that the values reported in Table 2 for static and fatigue onset were at 
least congruent with those data. 
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CONCLUDING REMARKS 


The purpose of this work was to study damage onset, accumulation and its 
consequences on glass/epoxy multidirectional laminates using simple stress 
analysis tools and fracture mechanics. The relatively lower stiffness of glass 
fibers compared to graphite fibers and the higher strains involved in damage and 
failure processes have highlighted the considerable importance of matrix cracks 
parallel to the fibers. Glass-fiber laminates are useful in studying matrix 
cracking but they show incongruences when other failure modes, such as edge 
delamination, are superimposed. 

A simple way to model stiffness reduction due to cracked 90° ply matrix was 
verified and shown to apply over a wide range of crack densities using both 
F.E.M. analysis and experiment. 

A closed form expression for strain energy release rate associated with 
matrix cracking has been proposed and, following experiments, strain energy 
release rate resistance curves (R-curves) for two different layups of the same 
material have been shown to coincide fairly well. The strain energy release rate 
approach used in this study can help explain the observed saturation crack 
density for matrix cracking. 

Edge delamination was difficult to analyze using the conventional tools 
developed for graphite reinforced composites. The superposition of matrix 
cracking leads to an overestimation of the interlaminar fracture toughness. 
Further investigation in this field is needed. 
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APPENDIX A : Finite element analysis of matrix cracking 


Matrix cracking was modeled using the well known NASTRAN finite element 
computer program, both to obtain a comparison of eq.(1) with a numerical stress 
analysis, and to evaluate the use of a general purpose structural program in 
fracture analysis problems. 

The shaded area of fig. 4 was modeled for different crack spacings. The only 

way to model out of plane fibers with the NASTRAN program was to use three- 

dimensional elements. A single thin layer was modeled using 8-noded 

isoparametric elements with constraints representing a plane-strain situation. 

The size of the smallest elements near the crack tip was 1/18 of the ply 

thickness. In total 208 nodes were needed for the [0/90] laminate model. 

For the [0/90] laminate the prediction of laminate plate theory was exactly 
s 

matched for the case with no cracks while a small edge effect was observed in 

the 45° and -45° plies of the [±45/0/90] laminate. The result of this edge 

s 

effect on laminate stiffness was eliminated by normalizing each stiffness value 
obtained for different crack spacings with the corresponding value obtained for 
a model of the same length without cracks. 

In conclusion it has been demonstrated that a general purpose program like 
NASTRAN is suitable to model this problem, although the use of three-dimensional 
elements makes it rather time consuming and expensive. 
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APPENDIX B : Exact solution for longitudinal modulus 


The longitudinal modulus of a [0/90] laminate may be calculated from 

s 

classical laminate plate theory as 


and due to the symmetry of the extensional stiffness matrix [A] of a [0/90] 

s 

laminate 


E 


x 




(B2) 


Introducing the reduced stiffness matrix components of the material Q llt Q 22 , 
Q 12 we obtain 


- 1 c(q„ ♦ *q„) ( 03 ) 

where x, the independent variable for this problem, is the ratio E 2 /(E 2 ) 0 of the 
transverse Young's modulus in the cracked 90° ply over its original value. 

If x=1 , eq. B3 gives E=E 0 , the stiffness of the undamaged laminate, while for 
x=0 it gives E=E^ used in eq. 3* 

Introducing the elastic properties of T190/x751/50 glass/epoxy, the maximum 
difference between eq. B3 and eq. 3 is less than 0.2 % and can be neglected. 
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Table I. Multiple regression parameters for fatigue damaging in a [±45/0/90j 

3 

glass/epoxy laminate. 

A. FATIGUE TESTS 


K a o 3l ^90 ~ az ^45 aj ^-45 a “*^b^ 

h ° y 4b 45 b ±45 b 0/90 


damage type parameter 

regression 

std. dev. 

undamaged condition 

a o 

.9945 

C\J 

•=T 

o 

o 

• 

90° ply matrix cracking 

a i 

.0674 

.0048 

+49° ply matrix cracking 

a 2 

.0415 

.0175 

-45° ply matrix cracking 

a 3 

.0238 

.0076 

+45/-45 interface delamination 

a H 

.321 1 

.0302 

0/90 interface delamination 

a s 

.1488 

.0357 

Coefficient of determination 

= 0.89 



B. 

. STATIC 

TESTS 


E 




Eo = 

a o " a i 

°90 b^ 


damage type parameter 

regression 

std. dev 

undamaged condition 

a 0 

1 .00 

.004 

90° ply matrix cracking 

a i 

.11 

.005 

0/90 interface delamination 

a 5 

.10 

.039 

Coefficient of determination 

= 0.92 
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TABLE II. Measured critical strain energy release rate for T1 90/x751 /50. 


damage 


load condition 

Gc 

90° ply matrix 

cracking 

static tension 

120 

J/m 

90° ply matrix 

cracking 

tension-tension fatigue 

30 

J/m : 

0/90 interface 

delamination 

static tension 

460 

J/ni 

0/90 interface 

delamination 

tension-tension fatigue 

110 

J/m : 
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EDGE DELAMINATION 



MAT«X CRACKS 
/ 45° PLY 


LOCAL DELAMINATIONS 


SPECIMEN EDGE 


Fig.1 : Photograph of a quasi isotropic [±45/0/90] glass-epoxy specimen, during 
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Fig. 2 : Schematic picture of possible damage modes of quasi isotropic layup. 
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Fig. 3 


Damage accumulation in a typical fatigue test on a [±45/0/90] specimen. 

s 

R=0 . 1 , o =160 MPa, Cycles to failure are 792300. 
max 
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Fig. 5 : Stiffness loss in static and fatigue tests on [0 /90 ] glass/epoxy 

m ns 


laminate. 






CRACK DENSITY*CRACK WIDTH 90 3 PLY EFFECTIVE TRANSVERSE MOOULUS . 

Fig. 6(a) Fig. 6(b) 


if) 



Fig. 6(c) 


Fig. 6 : Evaluation of 90° ply matrix cracking effect on generic laminate 
stiffness. 
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EQUATION (3) 



[+45/-45/0/90] c LAMINATE 




mnu 



CRACK DENSITY # CRACK WIDTH: D#w [No . #mm # [mm] 


Strain Energy Release Rate versus matrix crack density for different 


laminates . 





Comparison of Strain Energy Release Rate from analytical model (eq.8) 
and finite element analysis in a [0/90] glass-epoxy laminate. 
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